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ABSTRACT 

We use Minkowski Functionals (MF) to constrain a primordial non-Gaussian contribu- 
tion to the CMB intensity field as observed in the 150 GHz and 145 GHz BOOMERanG 
maps from the 1998 and 2003 flights, respectively, performing for the first time a joint 
analysis of the two datasets. A perturbative expansion of the MF formulae in t he limit 
of a w eakly non-Gaussian field yields analytical formulae, derived by Hikage et al.l 
(|2006D . which can be used to constrain the coupling parameter /nl without the need 
for non-Gaussian simulations. We find —1020 < /nl < 390 a t 95% CL, significantly 
improving the previous constraints by iDe Troia et al.l ()2007[ ) on the BOOMERanG 
2003 dataset. These are the best /nl limits to date for suborbital probes. 

Key words: Cosmology: early Universe - cosmic microwave background - methods: 
statistical - analytical 



2 P. Natoli et al. 



1 INTRODUCTION 

Detection of non-Gaussian signals in the Cosmic Microwave 
Background anisotropy pattern can be of significant help 
in discriminating between different inflationary models. The 
most simple inflationary models based on single-rolling 
scalar fields predict very small deviations from Gaussian- 
ity that cannot b e usefully constraine d by present day ex- 
perimental efforts l|Bartolo et al.|[2004h . However, multi-field 
infiationary models and other alternative scenarios allow 
for more relevant non-Gaussian contribution that could be 
in principle detected by current and forthcom ing missions 
(jBernardeau fc UzanI I2OO2I : iLvth et all l2003l l. In this pa- 
per we consider only the so-called local form for primor- 
dial non-Gaussianity, which can be parametrized by means 
of a quadratic terrn in Bardeen's curva l ^ure perturbations 
$ (ISalopek fc Bondl [l990l: iGangui et all Il994l : IVerde et all 
l2000l : lKomatsu fc Spergelll200ll ): 



<I>(x) = <I>g(x) + /nl [$g(x)' - ($g(x) 



(1) 



where "1>g is a zero mean, Gaussian random field and /nl 
is the coupling parameter that characterizes the ampli- 
tude of primordial non-Gaussianity. At present, the most 
stringent limits on /nl are derived from the WMAP five 
year analysis at —4 < /nl < 80 (95% CL) using an 
opti mal (i.e. minimum -variance) bispectrum based estima- 
tor (jSmith et al.ll2009l ). Many teams have further analysed 
the WMAP dataset to yield constraints on /nl using a 
plethora of t e sts, in c luding wavelet based e s timators: see e.g. 
ICurto et al.l (|2009l ): iRudiord et all (|2009l ): IPietrobon eFal] 
(|2009t ) and references therein. All /nl hmits to date are 
compa tible with a Gaussian hypothesis. lYadav fc Wandeltl 
(2008) claimed a measure of a positive /nl at above 99.5% 
CL in the WMAP three- year data using a bispectrum based 
statistics; however, their claimed signal ha s not been con- 
flrmed by the WMAP five-year analysis ()Komatsu et al.l 
l2008l : ISmith et al.ll2009l ). 

On the other hand, several groups have also investi- 
gated specific signatures in the WMAP data, typically in- 
duced by low resolution features such as anomalous spots, 
reporting high significanc e yet u nmodcUed detection of 
non-Gaussia.nity C Crcmin cUi et al.l [200 7: Cruz ct al._ 2007 ; 
Vielva et all l2004r lEriksen et al.l [2OO7I : iPark et al.1 120071 : 
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20091 ). 



CMB suborbital experiments hav e also delivered /n l 
constr aints, p articularly MAXIMA llSantos et al. 2003h. 
VSA (|Smith et al. 20^ ), BOOMERanG (iDe Troia etahl 
I2OO7I ) and ARCHEOPS jcurto et al.ll2008h . Although such 
limits are weaker than those based on WMAP, they probe 
a range of angular scale that will not be accessible to space 
borne observation until the advent of Planck Q. Among 
suborbital probes, iDe Troia et al.l (|2007l ) set the most strin- 
gent /nl constraints to date at —800 < /nl < 1050 (95% 
CL) from BOOMERanG 2003 (hereafter BOS) dataset using 
a pixel space statistics based on Minkowski functionals. Such 
constraints were obtained using a reference Monte Carlo set 
composed of non-Gaussian CMB maps. 

In this paper we revisit the /nl analysis of the 
BOOMERanG dataset. We employ a larger dataset that 
also includes the BOOMERanG 1998 (hereafter B98) data, 

^ http:/ /www.rssd.esa.int/index.php?project=planck 



allowing for a larger sky coverage and improved signal to 
noise. Furthermore, we apply a different, harmonic based, 
Minkowski functional code that overcomes a weakness of 
the previous BOS analysis, which used a flat sky approxi- 
mation to compute the functionals. Finally, we employ the 
perturbative approach developed bv iHikage et al.l (|2006l ) to 
qua ntify the contr i butio n of primordial non-Gaussianity to 
MF. iHikage et al.l (|2008!) successfully applied the perturba- 
tive method to WMAP data without need of a large set of 
non- Gaussian simulations. 

The plan of this paper is as follows. In Section [2] we 
briefly describe the BOOMERanG experiment and the two 
datasets it has produced as well as our data analysis pipeline. 
In Section[3]we apply the perturbative formulae to compute 
the MFs of the data and Gaussian Monte Carlo simulation 
maps. Furthermore in Section U we constrain /nl and in 
section [5] we draw our main conclusions. 



2 THE B98 AND BOS DATASETS 

BOOMERanG was launched for the first time from Antarc- 
tica in December 1998. It has observed the sky for about 
10 days, centering a target region at RA ~ 5/i and DEC 
~ —45° that is free of contamination by thermal emission 
from interstellar dust. BOOMERanG mapped this region 
scanning the telescope through 60° at fixed elevation and 
at constant speed. At intervals of a few hours the telescope 
eleva tion was change d in order to increase the sky cover- 
age ICriU et al.ll2003l ). The survey region aimed at CMB 
observations is ~ 3% of the sky or ~ 1200 square degrees 
(|Ruhl et al.ll2003l ). The data were obtained using 16 spider- 
web bolometric detectors sensitive to four frequency bands 
centered at 90, 150, 240 and 410 GHz. Here we restrict our- 
selves to the 150 GHz data that have the most advanta- 
geous combination of sensitivity and angular resolution to 
target the CMB fluctuations. The analysis of B98 dataset 
produced the first high signal to noise CMB maps at sub- 
degree resolution and one of the first high confidence mea- 
surements of the fi rst acoustic peak in the CMB anisotropy 
angular spectrum l|de Bernardis et al.ll20od ). The Gaussian- 
ity of this datas et has been constrained in both pix el and 
harmonic space jPolenta et al.ll2002l : iDe Troia et aL| [2003). 

The BOS experiment has been flown from Antarc- 
tica in 2003. Contrarily to B98, BOS was capable of 
measuring linear polar i zation other than to t al intensity 
(IPiacentini et al.' '20061; 'Mo ntrov ei~^ I2OO6I : Ijones et all 
2006; MacTavish et al. 2006). It has observed the microwave 
sky for ~ 7 days in three frequency bands, centered at 145, 
245 and 345 GHz. Here we use only the 145 GHz data, for 
reasons analogous to B98. These have been gathered with 
polarization sensitive bolometers (PSB), i.e. bolometers sen- 
sitive to total intensity and a combinatio n of the two Stoke s 
linear polarization parameters Q and U (| Jones et al.ll2003 ). 
The analysis of th e dataset has produced high quality maps 
(|Masi et al.ll2006l ) of the southern sky that have been con- 
veniently divided in three regions: a "deep" (in terms of 
integration time per pixel) survey region (~ 90 square de- 
grees) and a "shallow" survey region (~ 750 square degrees), 
both at high Galactic latitudes, as well as a region of ~ 300 
square degrees across the Galactic plane. The deep region is 
completely embedded in the shallow region. 
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In this paper we apply a pixel mask to select a larger 
effecti ve sky region than the one used in iDe Troia et al.l 
We have been extremely careful in choosing this 
sky cut, rejecting regions potentially contaminated by fore- 
ground emission, which shows up clearly in the B98 higher 
frequency maps, and pixels falling too close to the edge of the 
survey region, which exhibit low signal to noise and potential 
visual artefacts. The final cut we use covers about 980 square 
degrees or 2.4% of the sky. This shoul d be compared with th e 
~ 700 (1.7% of the sky) employed for lDe Troia et al.1 l|2007l ). 
which only used BOS, and with the 1.2% and 1 .8% of the sky 
selected re spectively for the B98 analyses of IPolenta et al.l 
l|2002l ') and lOe Troia et al.1 (|2003h . The sky mask employed 
in this paper is the largest ever used for BOOMERanG non- 
Gaussianity studies. 

We analysed the temperature ( T) data m ap reduced 
jointly from eight PSB at 145 GHz l|Masi et a l. 2006) for 
the BOS dataset and the T map obtained from the best 
three of the six 150 GHz channels for B98. While we do 
not consider here the Stokes Q and U polarization maps, 
the BOS temperature map has been marginalized with re- 
spect to linear polarization. The maps have been produced 
with ROMA, an iterative generalized least square (GLS) 
map making code specifically tuned for BOOMER anG anal- 
ysis (jNatoh et al.ll200ll : lOe Gasperis et al.ll2005h . To work, 
the GLS map maker needs to know the detectors' noise 
power spectral densities, which is estimated directly from 
flight data using an iterative procedure. In the case of BOS, 
where cross-talks among PSB located in the same feed horn 
are significant, we have also estimat ed the corresponding 
noise cross-spectra ()Masi et al.ll200^ ). The timelines have 
been carefully fiagged to exclude unwanted data; for B98, 
only the more conservative part of the scan surveye d at 1 
degree per second (d.p.s.) is used (|Crill et al.1 |200S| ) while 
for both datasets we have flagged all of the turn-around 
data. Once the B98 and BOS maps are produced, we obtain 
a single data map by noise weighting the two. In doing so 
we treat the residual noise left in the map as white. This 
choice is motivated by a property of the GLS map making 
procedure, which is very effective in suppressing the level 
of noise correlations in the data. The noise level in the B98 
map roughly equals that in the BOS shallow region: at 6.7' 
the noise r.m.s. is about 40 fiK per pixel. (While the B98 
flight actually lasted longer than BOS, we consider only three 
channels and the 1 d.p.s. part of the scan here.) The noise 
r.m.s. in the BOS deep region is ~ lO/xK for 6.7' pixels. The 
joint B0S/B98 map we obtain is shown in Fig. [1] in the sky 
cut employed for the analysis hereafter. 

To probe CMB non-Gaussianity it is important to 
keep under control contaminations from astrophysical fore- 
grounds, whose pattern is markedly non-Gaussian. In the re- 
gion selected here, foreground intensity is know n to be neg- 
ligibl e with respect to the cosmological signal l|Masi et al.l 
I2OO6I ) . We have masked all detectable sources in the observed 
fleld. To assess the robustness of our tests of Gaussianity we 
used a set of 1000 Monte Carlo simulation maps that mimic 
both the BOS and the B98 data. To produce these simula- 
tions, the following scheme is employed. The Gaussian CMB 
sky signal is simulated using the cosmo logical parameters es - 
timated from the WMAP 1-year data (jHinshaw et al.ll2003l l 
which flts well the BOOMERanG temperature power spec- 
trum. This signal is smoothed according to the measured 



optical beam and synth esized into a pixeli zed sky map, us- 
ing HEALPix routines (|G6rski et aLlbOOSl ). Taking into ac- 
count the BOS and B98 scanning strategy, the signal map 
is projected onto eight timestreams, one for each 145 GHz 
detector, for BOS and onto three timestreams for the B98 
150 GHz channels we consider here. Noise only timestreams 
are also produced, as Gaussian realizations of each detec- 
tor's noise power spectral density, iteratively estimated from 
flight data as explained above, taking into account cross 
talks between detectors in the case of BOS. The signal and 
noise timelines are then added together. To reduce the sim- 
ulated timelines, we follow the same steps performed when 
analysing real data: the timelines are then reduced with the 
ROMA mapmaking code replicating the actual flight point- 
ing and transient flagging, to produce T maps jointly from 
three B98 channels and T, Q and U maps jointly from all 
eight BOS channels. We enforce that the map making proce- 
dure is applied to simulation and observational data follow- 
ing the same steps. 

It is worth noticing that in this paper the B98 and BOS 
data have been used to produce a joint map for the first 
time. 



PERTURBATIVE APPROACH TO 
MINKOWSKI FUNCTIONAL FOR A 
WEAKLY NON-GAUSSIAN CMB FIELD 



In the previous paper (|De Troia et al.ir2007l ) we have also ap- 
plied simple pixel beised analysis (specifically, the normalized 
skewness and kurtosis) to the BOS observed field. Here we 
restrict ourselves to three Minkowski Functionals generally 
defined in two-dimensional maps: fraction of area Vb, total 
circumference Vi, and Euler Characteristic V2. We measure 
the MFs for CMB temperature maps above the threshold 
density v, defined as the temperature fluctuation / = AT/T 
normalized by its sta ndard devia ti on a = (f ^)^^^. Follow- 
ing th e formalism by iMatsubaral (|2003l ) and iHikage et al.l 
(|2006l ). we can write the analytical formula for the fc-th 
Minkowski functional of weakly non-Gaussian fields as 



UJ2 



k{k~ 1) 
6" 



V2 



(2) 



(TO + 0((To) 



where Hn{v) are the Hermite polynomials and ujk ~ 
7r'''''^r(fc/2 + 1) gives uo = 1, tJi = 2, and 0^2 = tt. The 
S^''\i = 0, 1, 2) are skewness parameters, defined by 



q(0) _ if) 



3 (,f (VV)) 



5(2) ^ 3 ((V/)-(V/)(V^/)) ^ 



The variances o-j (j — 0, 1) are calculated from Ct as 



(3) 
(4) 
(5) 

(6) 
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Figure 1. The CMB field as seen in the B2K + B98 map, in th e sky cut used for the analysis presented here. The solid line shows the 
boundary of the region taken in consideration for the analysis in lDe Troia et aL 1 1I2OO7I') . 



where Wt is a window function that includes the experi- 
ment's effective optical transfer function (assumed circularly 
symmetric) and low-^ and high-^ cut-off as well as the filter 
function due to pixelization effects. Expanding the skewness 
parameters into spherical harmonics a nd using the reduced 
bispe ctrum bi^e^i^ as a function of /nl (|Komatsu fc Spergell 
[200J), we get: 



-AO) 



?(2) 



27r(T, 



E 



'ie2i3btie2e3We,We2We^, 
[^1(^1 + (^2 + 1) 



3 



(7) 

(8) 



{[£1(4 + 1) +4(4 + 1) 



-4(4 + 1)]4(4 + 1) + (eye.)} 

X Ii^e2i3 bi^i2i3 We^ Wi2 Wt^ , 



where 



^hhh 



(2Zi -1-1) (2^2 + 1) (2^3 + 1) 



47r 



I2 




(9) 
(10) 



(11) 



In the above theoretical predictions we assume a ACDM 
model with the cosmological parameters at the maximum 
likelihood peak from WMAP 1-year data (jSpergel et al.l 
L2OO3): = 0.043, f^cdm = 0.21, Oa = 0.74, Hq = 72km s'^ 
Mpc~^, Us — 0.96, and r = 0.11. The amplitude of primor- 
dial fluctuations has been normalized to the first peak ampli- 
tude of the temperatu re power spectrum , £{£+ l)Ci/{27v) = 
{7A.7fiKf at £ = 220 (|Page et al.ll20o3 ). 



We compute the MPs of the pixelized maps by integrat- 
ing a combination of first and second angular derivatives of 
th e temperature over the sky, as described in Appendix A.l. 
of lHikage et~al] l|2006i 'l. The threshold density v is set in the 
range —3.6 to 3.6, assuming 19 evenly sp aced grid points. For 
our analysis we use maps at HEALPix l|G6rski et al.ll2005h 
resolution of ~ 13' (iVsido = 256) and ~ 7' (iVnsidc = 512). 



4 CONSTRAINTS ON PRIMORDIAL 
NON-GAUSSIANITY 

We define a "joint" estimator by grouping the 14 's in a sin- 
gle, 57 elements data vector Vj = {Vo{i' — — 3.6), Vo(!^ — 
-3.2),..., Vb(i/ = 3.6),Vi{u = -3.6),...,V2{iy = 3.6)}. We 
want now to constrain the /nl parameter and estimate its 
best fit value. Starting from analytical formulae we can cal- 
culate the non- Gaussian part of the MPs using the equation 
©, i.e.: 

A1/j(/nl) = Vj(/nl) - Vj{fyi. = 0). (12) 

We can then estimate our final non-Gaussian predictions as: 

V'j(/nl) = 'K7(/nl = 0) + A1/7(/nl), (13) 

where Vj(/nl = 0) is the average MP computed from 
our Gaussian Monte Carlo simulations. The reason for this 
choice is that the Monte Carlo average provides an accu- 
rate estimate of the MPs, accounting for instrumental and 
coverage effect. Pinally we perform a chi-square analysis by 
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Figure 2. The plots show, from top to bottom, Vq, V\ and V2 computed from the BOOMERanG data (filled circles), from analytical 
formulae (solid line, here computed with the best fit /nl value for each functional) and from Gaussian Monte Carlo simulations (dashed 
lines). The columns refer to different choices of the low £ cut (see text), £niin = 6, 12 and 20 (from left to right). The error-bars for the 
data are derived as the 1-cr standard deviation of our Monte Carlo simulations. The pixel size of the map is ~ 7' (HEALPix Afgijje = 512). 



measuring 



ulations. We have verified that, when computing its ma- 
trix elements, one needs to take into account the correla- 
Bgs+Boa^^-j _ Vj{f-M-L)]C J j, [V^'^*^^"'' — Vj'(/NL)],tions among different functionals not to incur in biased con- 
straints. In Fig. [2] we plot each MF of the B98 and BOS data 
(14)compared with the theoretical predictions with the best fit 
value of /nl for each MF. The error-bars are derived as l-u 
deviations estimated from 1000 Gaussian maps with corre- 
lated noise. 



where vf^^+^°^ denote the MFs for the joint B98 and BOS 
map. This expression is used to derive constraints for /nl 
and for our goodness of fit analysis. The full covariance ma- 
trix Cj^ji is also estimated from Gaussian Monte Carlo sim- 



We study the effect that neglecting the contribution of 
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a range of multipoles I has on this analysis. A low £ cut is 
necessary since we are dealing with data from a suborbital 
experiment, which has not been designed to measure large 
angular scales. These cannot be constrained properly, firstly 
because of the limited angular extension of the region sur- 
veyed, and secondly because timeline filtering is applied to 
the data to suppress contribution from low frequency noise 
and systematics. The filters are applied during th e map mak- 
ing stage at ~ 70 mHz both for BOS and for B98 l|Masi et all 
I2OO6I : ICriU et al.ll2003l ). While we apply the same fihers in 
our simulations, the amount of low i power in the latter is 
somewhat different from those exhibited by the data. This 
happens because the GLS map maker is more efficient in 
recovering the large scale structure from the simulations, 
where we only add Gaussian noise with random phases, than 
from the data where the noise has a more complex structure. 
To account for this effect, we exploit one degree of free- 
dom allowed by the harmonic analysis to MF here pursued: 
specifically, we set to zero all the power in the map below 
(■-mi-n- The left to right panels of Fig.[2]refer to three different 
choices: £min =6, 12 and 20 respectively. The imin = 6 cut is 
the natural one that would arise due to limited sky cover- 
age but in this case the MFs from data and simulations do 
not agree well for Vi and V2 (see Fig [2]) . The agreement is 
much better for fmin = 12 and Imiu = 20 with no appreciable 
difference between the two. This comes with little surprise, 
since a telescope scan speed in the range 0.5° /s to l°/s, both 
of which have been employed in the dataset we consider, ef- 
fectively high pass filters the data in the range 10 < £ < 20. 
(Note, however, that timeline filtering has an anisotropic ef- 
fect on the sky, due to the nature of the scanning strategy 
employed for BOOMERanG.) 

On the other hand, it is advantageous to consider also a 
high i cut ^max to probc how the decreasing signal to noise 
level can affect /nl constraints. For the data, this can be 
done by varying the HEALPix resolution parameter A'sidc 
which is linked to Iraax in the spherical harmonic transform. 
We focus in the following on both A'aide = 256 with fmax = 
512 and A'side = 512 with ^max = 1000. Our dataset is signal 
dominated a t ^ ~ 500 while b egins to be noise dominated 
at f ~ 1000 (| Jones et al.ll2006l l. 

In Fig. [3] we show the analytical non-Gaussian correc- 
tions AVj (eq. [12]) for each MF compared to the residuals 
obtained by subtracting to the B98/B03 data MFs their 
Monte Carlo average, that is, V"B98+bo3 _ Vj(^f^^ ^ q). The 
error-bars represent the la error estimated from 1000 Gaus- 
sian Monte Carlo simulations. The analytical residuals are 
computed using the best fit value of /nl as obtained by 
minimizing the in the equation (|14|l . albeit this is done 
separately for each MF, ignoring (only for the sake of this 
plot) correlations among different functionals. The analyti- 
cal AVj in Fig. [3] are normalized to the maximum of their 
Gaussian part, while the data points are normalized to the 
maximum of the Monte Carlo average. We show results both 
for A'aide = 256 (left) and for A^sidc = 512 (right) and for a 
low-multiple cut at ^min = 6, 12 and 20 (left to right). The 
agreement of the residual plots suggests that it is safer to 
adopt the most conservative cut at Imin = 20. Since increas- 
ing further the cut does not appear to yield further advan- 
tage, we focus on a conservative choice ^min = 20 for our 
final analysis. 

Table [1] lists the confidence intervals for /nl estimated 



from the BOOMERanG data at different lyain and ^max 
thresholds. The results are obtained taking into account the 
full covariance matrix of the Vfc's, as expressed in the equa- 
tion (O above. The MFs computed at different HEALPix 
resolution encode different statistical information on the un- 
derlying field and thereby combining two sets of MFs im- 
prove the limits on /nl. We build a global covariance matrix 
to take also into account correlations among the two sets. 
We repeat this process for each value considered. 

In the conservative case of fmin = 20 using the "com- 
bined" estimator, our analysis yields —670 < /nl < 30 at 
1-cr level and -1020 < /nl < 390 at 95% CL, while the min- 
imum (best fit) value of /nl is at —320 . The previous analy - 
sis performed only on the B03 dataset (|De Troia et al.|[2007f ) 
produced limits weaker than those obtained in this paper by 
a factor ~ 1.4. The improvement is due to a combination of 
both enhanced sky coverage and the use of the analytical 
estimator. The B98 dataset alone has never been subject 
so far to /nl analysis. We have checked that this dataset 
too is compatible with /nl = at l-cr with A/nl = 470 
(A/nl = 930 at 2-cr). 

We can also quantify the cost of imposing a low I cut 
to the data. In fact, had we not considered an effective Imin 
value, one would expect to reduce the confidence interval on 
/nl by ~ 1.6 (c.f. the last row in Table[T] obviously obtained 
not from the data but from a simulation containing a low 
resolution pattern). In practice, this could be obtained by 
adding to the dataset a low resolution CMB field coming, 
e.g., from the WMAP data. While this would give us tighter 
constraints on /nl , we prefer to focus here on the limits one 
can derive from the BOOMERanG data alone. Note also, 
that a diminished sensitivity to low resolution features is 
a characteristic common to most -if not all- of the subor- 
bital experiments. The accurate measurement of the CMB 
at low and high multipoles with one single experiment is 
rather a prerogative of space born missions, which enjoy the 
necessary stability and long term integration capability. Our 
analysis is the first (to our knowledge) to explicitly take into 
account this effect for a suborbital experiment. To explain 
the significant broadening of /nl constraints caused by the 
low £ cut, one can note that for the underlying ("local") 
form of non-Gaussianity we are probing here, the low multi- 
poles are actually very important. In fact, most of the signal 
in the reduced bispectrum be-^e^e^ lies in "squeezed" £-space 
triangles, with one side much smaller than the other two. 
When probing non-Gaussianity one is basically comparing 
signal at the lowest multipole with two of the highest mul- 
tipoles. As a result, S/N increases as ^max/^min so one can 
either increase ^max for a given Imin (which explains, e.g., 
the improvement of WMAP over COBE and the forecasted 
improvement of Planck over WMAP), or reducing Imin for 
a given f^ax. 



5 CONCLUSIONS 

We have analysed data from the BOOMERanG experiment, 
combining for the first time the temperature maps of the 
1998 and 2003 campaigns, to constrain a non-Gaussian pri- 
mordial component in the observed CMB field. We focused 
on Minkowski Functionals, comparing the data to analytical 
perturbative corrections in order to get constraints on the 
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Figure 3. Comparison of the three MFs residuals for B98/B03 temperature data (filled circles) to the analytical predictions with the 
best fit value of /nl for each functional (solid lines). The analytical predictions are normalized to the maximum value of the Gaussian 
part while the data points are normalized to the maximum of the Monte Carlo average. The left figure is for HEALPix Afgi^e = 256 
resolution while the right figure is for N^idc = 512. From left to right in each figure, we show the l^nin = 6, 12 and 20 cases. The error-bars 
represent the standard deviation at l-cr estimated from 1000 Gaussian Monte Carlo simulations. 
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Table 1. The confidence intervals for /nl estimated for the 
B03/B98 data with different low and high £ cut values (see text). 
We show both the la and 2a confidence interval. The last two 
columns are obtained from a combined analysis of data for the 
two £max values. The last row is derived using a simulation with 
no effective low £ cut and shows the improvement that could 
be obtained if a low resolution pattern had been present in the 
BOOMERanG field. 



non- linear coupling parameter /nl- We have used a set of 
highly realistic simulation maps of the observed field gen- 
erated assuming a Gaussian CMB sky, since the formal- 
ism we have adopted does not require non-Gaussian sim- 
ulation maps. We studied the effect that the lack of low 
resolution CMB features in the BOOMERanG data has on 
/nl constraints. We find -670 < /nl < 30 at 68% CL and 
-1020 < /nl < 390 at 95% CL. These limits are signifi- 
cantly better than those published in a previous analysis lim- 
ited to the BOOOMERanG 2003 data (-800 < /nl < 1050 
at 95% CL), represent the best results to date for suborbital 
experiments and probe angular scales smaller than those ac- 
cessible to the WMAP. 
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